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ABSTRACT 

We have selected a sample of X-ray-emitting active galactic nuclei ( AGNs) in low-mass host galaxies 
(~5xl0^- 2x10^" Mq) out to z~l. By comparing to AGNs in more massive hosts, we have found 
that the AGN spatial number density and the fraction of galaxies hosting AGNs depends strongly on 
the host mass, with the AGN host mass function peaking at intermediate mass and with the AGN 
fraction increasing with host mass. AGNs in low-mass hosts show strong cosmic evolution in comoving 
number density, the fraction of such galaxies hosting active nuclei and the comoving X-ray energy 
density. The integrated X-ray luminosity function is used to estimate the amount of the accreted 
black hole mass in these AGNs and places a strong lower limit of 12% to the fraction of local low-mass 
galaxies hosting black holes, although a more likely value is probably much higher (>50%) once the 
heavily obscured objects missed in current X-ray surveys are accounted for. 
Subject headings: galaxies: nuclei - galaxies: active - X-rays: galaxies 
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1. INTRODUCTION 

Active galactic nuclei (AGNs), the manifestation 
of accretion onto massive black holes (MBHs), have 
been recognized as a critical ingredient in galaxy for- 
mation and evolution. The demography of local 
galaxies suggests that most - perhaps all - massive 
galaxies host MBHs at their centers and that MBH 
masses are correlated with the ga l axy bulge prop- 
erties (iKormendv fc Richstond 119951: iMagorrian et al.1 



19981: iGebhardt et all 120001 : iFerrarese &: MerrittI 120001 : 



Haring fc Rixl l2004l ). implying the coevolution of the 
galaxy and MBH. The good match between the 
local BH mass density and the mass density of 
AGN relics further suggests that all massive galaxies 
have exp erienced an AGN phase during their evolu- 
tion fe.g. lAller fc Richstond 120021: IShankar et all l2004l : 
iMarconi et al.1 I2004D . Energv feedback from AGNs to 
their host galaxies is invoked to explain different as- 
pects of massive galaxy evolution. AGNs may serve as 
the heating sources for co oling flows in clusters (e.g. 
iMcNamara fc NulsenI 120071 ) . They may suppress star 
formation in their host galaxies and cause them to mi- 
grate from the blue cloud to the red-sequence in the 
color - magnitude plot (ICroton et all 120061 : iNandra et al.1 
120071 : iGeorgakakis et al.ll2008l ) and t hey may account fo r 
"down-sizing" galaxy evolution (e.g. ICowie et al.lll996f ). 
As the most abundant population in the universe, low- 
mass (defined as stellar mass M* < 2x10^" Mq through 
this paper) galaxies act as the building blocks of massive 
galaxies. However, the current understanding of MBHs' 
role in low-mass galaxy evolution is limited, leaving us 
with some basic questions: Is the existence of BHs in low- 
mass galaxies as common as it is in massive galaxies? Are 
there two types of low-mass galaxy populations (i.e. ones 
with BHs vs. ones without BHs)? How many low-mass 
galaxies experience an AGN phase? All of these ques- 
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tions are related to a more basic astrophysical problem: 
what is the black hole occupation function (BHOF; the 
fraction of galaxies hosting either active or quiet BHs) in 
low-mass galaxies? 

Searching for BHs in low-mass systems offers an unique 
opportunity to extend the MBH-cr correlation for mas- 
sive galaxies, which not only tests the universality of the 
relation but also is required to understand the origin 
of the relation. Although the MBH-cr relationship has 
been confir med for low-mass systems in the local universe 
(jBarth et al. 2005 ) , some works suggest that the MBH-cr 
relationship may be replaced at low m ass by a similar re- 
latio n involving com pact stellar n uclei (|Wehner fc Harris! 
120061 IFerrarese et"a l. 2006; Ross a et al.ll2006D . 

A better understanding of primordial MBH seed 
growth in the early universe may also benefit from the 
study of BHs in low-mass galaxies, as such work could 
reveal aspects of the accretion mode in a low gravita- 
tional potential and low metallicity environment and of 
the relative importance of mass growth through accre- 
tion and merging processes. For example, in low-mass 
galaxies, the impulsive kick from anisotropic gravita- 
tional emission during BH-BH merg ers is thought to be 
strong enough to ej ect central BHs (|Favata et al.ll200i : 
iMerritt et"aLl l2004f ). Understanding AGN activity in 
low-mass galaxies can further constrain MBH seed for- 
mation theories. In a cold dark matter universe, the 
primordial MBH seeds form as re mnants of Population 
HI stars (e.g. iMadau fc Reesll200ll) or through the direct 
collapse of pre-galactic gas disks (|Lodato fc NataraianI 
[200 6) . The efficiencies in different formation scenarios 
predict a wide range in the local BHOF in low-mass 
galaxies which can vary from zero to unity, while the 
prediction of the BHOF in massive galax ies is invariably 
unity as constrained by observations (|Volonteri et al.l 
I200I . 

The current study of BHs in low mass galaxies 
is mostly limited to low redshift {z < 0.3). De- 
spite a variety of investigations, it is still unclear how 
many local low-mass galaxies harbor BHs. The dy- 
namical searches for MBH s have confirmed th e ex- 
istence of a BH in M32 (|Verolme et al.l I2002D . but 
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not in M33 (iMerritt et all 120011: JGebhardt et aLll2001h 
or NGC205 (jValluri et atl l2005l ). Searching for ac- 
tively accreting BHs (i.e. AGNs) has provided a 
stronger lowerUmit and more complete view of local 
BHOFs as a function of the system mass through 
optical emission l ines (e.g. iFiliPDenko fc Sargent) Il989l: 
Maiolino fc Riekel[l 995: Ho et al.lll995HKauffmann et all 
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20071: [Shields et al 



2003at iGreene fc Hq ,2004 l2007cl: iDecarU et all l2007[ 



selections ( Satvapa l ct al. 2007, '200 8D and X-r ay emis- 
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sion fe.g. fGallo et al. 2007 : Ghosh etall 120081) . It has 
been shown that the fraction of galaxies hosting active 
BHs depends on both BH mass and galaxy stellar mass, 
peaking at an intermediate m ass range and falling to- 
ward higher and lowe r mass dK auffma nn et al.l l2003al : 
iHeckman ct al. 2004 : Greene fc Ho 2007bi r 

As the local BH masses are most likely accumulated 
at high redshift, the search for AGN activity in low- 
mass galaxies at high redshift should provide indepen- 
dent and possibly better constraints on the local BHOF. 
Such study is also significant for understanding the evolu- 
tion of nuclear activity in low-mass galaxies and the role 
of AGNs in low-mass galaxy evolution. With the advent 
of deep Chandra and XMM-Newton X-ray surveys, we 
have searched for X-ray emitting AGNs in low-mass host 
galaxies out to z ^ 1. In this paper, we first describe 
the identification of AGNs in low-mass galaxies (see § [2] 
and §[3]). We then present the l/Knax method to correct 
for incompleteness in § |4] and study their spatial number 
density and X-ray luminosity function in § [S] In § |6l we 
discuss the local BHOF in low-mass galaxies constrained 
by our study of high redshift AGNs in such galaxies. Our 
conclusions are presented in §[71 Throughout this paper, 
"low-mass" refers to normal galaxies or AGN host galax- 
ies with stellar mass M* < 2x10^° M© and "massive" 
indicates those with stellar mass Af, > 2x10^" M©. We 
adopt a cosmology with Ha—7Q km s^^ Mpc^^, il,n=0.3 
and f7A=0.7. All magnitudes are defined in the AB sys- 
tem. 

2. DATA: X-RAY SURVEY FIELDS 

We have searched for AGNs in low-mass host galax- 
ies in five Chandra and XMM-Newton fields, includ- 
ing the All- Wavelength Extended Groth Strip Interna- 
tional Survey (AEGIS), the Chandra Deep Field-North 
Survey (CDF-N), the Chandra Deep Field-South Survey 
(CDF-S), the Chandra Large- Area Synoptic X-Ray Sur- 
vey (CLASXS) and the XMM-Newton Large Scale Struc- 
ture Survey (XMMLSS). Table [T]hsts the properties of 
these five fields, including the area, the limiting hard X- 
ray flux, the available optical/near-IR photometry, the 
definition of secure (multiple- line) spectroscopic redshifts 
and the associated references. 

All the X-ray fields are fully covered by optical/near- 
IR photometry. For the CDF-N, CDF-S and CLASXS 
fields, spectroscopic observations have been obtained for 
X-ray sources, while the redshifts of X-ray objects in the 
remaining two fields are obtained by matching X-ray cat- 
alogs to galaxy redshift survey catalogs. Table [2] summa- 
rizes the total number of X-ray sources and of spectro- 
scopic targets. The search radii for optical counterparts 
to the X-ray sources are 2.0" and 2.5" for the Chandra 
and XMM-Newton fields, respectively. The optical coun- 
terparts are identified using _R-band catalogs, which gen- 



erally provide the deepest observations. If multiple opti- 
cal objects within a search aperture are present, the clos- 
est one is defined as the optical counterpart. The fraction 
of X-ray objects with multiple optical sources within a 
single aperture is only ^ 10%. Therefore, the assump- 
tion of adopting the closest one as the optical counterpart 
should not affect our conclusions. We have limited our 
study to objects with secure spectroscopic redshifts (see 
§ 131) and thus the majority of the optical counterparts 
are brighter than 24 in the _R-band. At R <24, the sur- 
face density of gala xies and stars is about 16.6 arcmin^^ 
(jCapak et al.ll2004 ). The probability for chance superpo- 
sition between X-ray and optical sources is 6% and 10% 
for the Chandra and XMM-Newton fields, respectively. 
Given a total of 32 low-mass AGN hosts in the Chandra 
fields and zero in the XMM-Newton field, we have esti- 
mated that only two objects are expected to have spuri- 
ous optical counterparts. Our sample only includes the 
X-ray sources with detected hard X-ray fluxes, defined 
in the energy range of 2-8 keV. The published 2-10 keV 
fluxes in the AEGIS and XMMLSS fields have been cor- 
rected to 2-8 keV on the assumption of a power-law pho- 
ton index of 1.0, which is the average of hard- X-ray se- 
lected objects based on the hardness ratio (jNandra et al.l 
[20051) . 

All optical type 1 AGNs are excluded as their nuclear 
radiation contaminates the host optical/near-IR light 
severely. Due to lack of access to the observed spectra, 
the definition of optical type 1 AGNs is not completely 
universal over all the fields. In the CLASXS field, type 
1 objects are defined by having an emission line FWHM 
> 1000 km s-\ while FWHM > 2000 km s^^ is adopted 
for the CDF-N and CDF-S fields. For the AEGIS and 
XMMLSS field, we have downloaded the spectra and 
classified type 1 objects using FWHM > 1000 km s~^. In 
the AEGIS field, some fraction of type 1 objects is still 
included in the sample as the DEEP2 spectral coverage 
misses the permitted lines (Ha, H/3 and MgII2800A) in 
certain redshift windows (~0.2-0.3 and ~0.6-1.2). Given 
a broad line AGN fraction of 16% within broad-line- 
detectable redshift ranges (0-0.2 and 0.3-0.6) and a to- 
tal of eight AGNs in low-mass hosts within broad-line- 
undetectable redshift ranges in the AEGIS field, only one 
object in our final sample may be misclassified as an op- 
tical type 2 AGN. 

3. SELECTION OF AGNS IN LOW-MASS HOST GALAXIES 

Objects are identified as active low-mass galaxies if 
they satisfy the following two criteria: 1.) stellar mass 
Af, < 2x10^" Mq, our definition of low-mass galaxies 
and 2.) hard X-ray luminosity i2-8kcV > iC^ erg s~^, 
indicating an active nucleus. 

We have measured stellar masses by comparing the ob- 
served SEDs (i.e., the photometric data included in Ta- 
ble [Tj) to 102168 s t ellar synthesis models produced by 
iBruzual fc Chariot! (|2003[ )'s code. As hsted in Table [3 
the stellar models span a wide range of parameter space, 
including metallicity, extinction, characteristic timescale 
of exponential star formation history, fraction of ejected 
gas being recycled and galaxy age. To account for the 
possible existence of low metallicity and young galaxies, 
we have included all six available metallicities and galaxy 
ages st arting at 10^ y rs. T he fit algorithm is similar to 
that of iBundv et al.l (|2006f ) who have used a Bayesian 
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Fig. 1. — Examples of the best-fit SEDs and the probability dis- 
tribution of the stellar mass of AGNs in low-mass galaxies in each 
field, where the solid line is the median stellar mass and two dotted 
lines indicate the 16 % and 84 % probability tails, respectively. 



technique as described in iKauffmann et all (|2003cD . In 
summary, for an individual object, the best-fit steUar 
mass is obtained for each model that corresponds to an 
age younger than the cosmic age at the redshift of the 
object. The associated x^ gives the probabihty (exp(- 
X^/2)) that this model represents the observed SED. 
The final probability distribution of the stellar mass is 
obtained by summing all probabilities within a certain 
mass bin. We adopted a bin width of 0.01 in logAf*, 
which on average contains 300 models. The median value 
of the stellar mass probability distribution is adopted as 
the final mass of a galaxy. Compared to the minimum x^ 
derived stellar mass, the mass obtained by this technique 
suffers much less from model degeneracies. The 68% un- 
certainty range of the derived stellar mass is defined by 
excluding the 16% tail at each end of the probability 
distribution. Note that this uncertainty mainly reflects 
the photometric errors. For the objects in the CDF-N 
and CLASXS fields, there are no published photomet- 
ric errors. We have adopted universally an uncertainty 
of 0.07 magnitude, which is roughly the sky noise for a 
7Tiij=24 object in these two fields. A small fraction of ob- 
jects in AEGIS, CDF-S and XMMLSS with very small 
photometric errors show large minimum reduced x^, in- 
dicating our stellar models are not able to produce the 
observed SEDs accurately. To estimate the uncertainty 
for these objects, we increased their photometric errors 
to 0.02 magnitude, which gives a reasonable minimum 
X^ (< 10) and larger uncertainty. 

We do not include other systematic errors for our 
measured stellar masses, such as the accuracy of the 
iBruzual fc Chariot! (j2003l) code itself and possible alter- 
native choices for the initial mass function. Our sample 
only includes narrow-line AGNs and thus AGN light con- 
tamination to the host emission should be small. For ob- 
scured A GNs, the scat tered nuclear light may be strong 
((Zakamska et al.l |2006[) . However, such scattered emis- 
sion may affect the host light importantly only in the 



> 10^^ erg s ^ AGN selection criterion is 



UV band. To avoid the dust emission from AGN dusty 
tori, we have not used photometry in the mid-IR. Fig. [T] 
shows examples of the observed SED superposed with 
the minimum-x^ stellar model and the probability dis- 
tribution of the stellar mass for each field. 

A mass cut of < 2x10^° Mq is adopted to define the 
sample of low-mass host galaxies. Selecting the sample at 
this mass threshold is of interest because the properties of 
galaxies and BHs may transition around this mass. Stud- 
ies of low-redshift galaxies have shown that galaxy prop- 
erties (star formation history, size and internal structure) 
show significant differences at the dividi ng stellar mass 
of 3x101° Mq (jKauffmann et al.ll2003bh . In addition, 
galaxies with stellar mass < ^lO"'^" M© often harbor com- 
pact stellar nuclei at their centers, which follow the MBH 
mass-bulge relati onships, but w hich arc rare in more 
massive galaxies (ICaroUo et al.l 11 998: Lainc et al. 200^ 
IWehner fc Harris! l2006t iFerrarese et al., ,2006 ). These 
stellar nuclei may be replacements for MBHs in low-mass 
systems where the gravitational potential is not suffi- 
ciently deep to form a BH. Furthermore, the fraction of 
galaxies hosting active BHs depends on the galaxy stel- 
lar mass jK auffma nn et al.ll2003at iHeckman et alll2004l : 
iGreene fc H o 20071^). In the most complete local AGN 
sample, the AGN fraction in galaxies fainter than Mb 
— -20 (^lO^^ MfTi) is on aver age half of the fraction in 
brighter ones (|Ho et al.lll997l ). 

The L2-8koV 

mainly based on the energy budget argument that star- 
forming galaxies rarely produce such high hard X-ray lu- 
minosities (e.g. Zezas et al. 1998). The resulting sample 
of AGNs in low-mass hosts contains 32 objects as listed 
in Table H) Each object is labeled by the field name 
followed by the sequence number in the X-ray catalog of 
each field (see Table[T]for references). Fig.[5]shows distri- 
butions of low-mass AGN properties, including redshift, 
host stellar mass, 2-8 keV rest-frame luminosity and X- 
r a.y to J?-band flux ratio. 

iBest et al] (|2005D have measured average ratios of BH 
masses to the galaxy stellar mass as a function of stel- 
lar masses for the SDSS galaxy and AGN samples (see 
their Fig.l), where the BH mass is measured through 
the stellar veloclity dispersion and the MBH-cr relation. 
By assuming that our X-ray selected AGN sample has 
the same BH-to-galaxy mass ratio as the SDSS AGN, 
we can estimate the Eddington ratio for galaxies with 
M*= 2x10^° Mq and L2-8keV=10'*^ erg s'^. By assum- 
ing the Lhoi/L 9_inw.v = 17(L2-iokev/10^^ erg s-i)0-43 
(jShankar et al.l [20(M) and L2-8kev/i2-iokcV=0.86 (as- 
suming a power law photon index of 1.0), we have an 
Eddington ratio of 0.007. As shown in §[4l although we 
can detect galaxies with stellar masses down to around 
10* Mq, our sample is complete down to stellar mass of 
10^'^ Mq out to redshift of 0.7. For this st ellar mass, 
the b lack hole mass is on average 3x10^ Mq (jBest et al.l 
|2005( ) and the corresponding Eddington ratio is 0.02 for 
i2-8kcV=10''^ erg s^i. As a comparison to those local 
low-mass AGNs, our sample may be not deep enough 
to include the local classical low-mass Seyfert galaxies 
similar to NGC 4395 or POX 52, but should include 
galaxies simil a r to some low -mass SDSS AGNs found by 
IGreene fc H"oI (12004'. '2007c]) . 

To better understand the role of AGNs in low-mass 
galaxy evolution, we will compare the sample of low-mass 
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Fig. 3. — The effective solid angle as a function of the hard X-ray 
flux for objects in AEGIS, CDF-N, CDF-S and CLASXS fields. 



Fig. 2. — Distributions of properties of AGNs in low-mass galax- 
ies: (a)— redshift distribution; (b) - distribution of host stellar 
masses where the dotted line indicates the completeness cut of 
10^-^ Mq; (c) - distributions of the rest-frame 2-8 keV X-ray lu- 
minosity; (d) - distribution of X-ray to i?-band fiux ratio. To be 
complete, our number density and X-ray luminosity function for 
AGNs in low-mass hosts are measured for those at 0.1 < z < 0.7, 
lO'''^ Mq < M. < 10i°-3 Mq and ij-kev > 1°*^ «^^g '^"^■ 
host AGNs to a sample of AGNs in massive host galaxies. 
The comparison sample is defined as stellar mass M^ > 
2x101'' Mq and La-skcV > lO'*^ 
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4. INCOMPLETENESS CORRECTIONS AND WEIGHTS 

4.1. The i/Knax method 

The spatial number density of AGNs as a function of 
AGN host mass and the X-ray luminosity function of 
AGNs in low- mass host s was calculated using the 1/T4iax 
method (|Schmidtlll968D . This calculation used those ob- 
jects with R < 24, L^!!*8kev > lO**^ erg s'^ and M, > 
5 X 10^ M© and in two redshift intervals of 0.1 < z < 0.4 
and 0.4 < z < 0.7. ^a-skeV i^ the rest-frame 2-8 keV flux 
assuming a photon index of 1.0 but not correcting for ex- 
tinction. The maximum volume over which an object is 
included in the sample is given by 



K, 



ll——az, 
dz 



(1) 



where [zjow , Zhigh] is the redshift range of interest and Q, 
is the solid angle covered by the X-ray survey at the flux 
level of the object. While ziow is always fixed to the low 
end of a redshift interval, the maximum redshift, Zhigh is 
defined as: 



7u- I - minfr'"^;!! limit In 
Zhigh — nmn^z^^jj^ , z^^ , Zp, 



limit .^limit\ 



(2) 

where Zj^J^^ is the high end of a redshift interval, z^™^ 
is the limiting redshift at which the observed X-ray flux 
reaches the limiting flux in a given field where the K- 
correction is determined by assuming a power law spec- 
trum with photon index of 1.0 and z^™'' is the limiting 
redshift where the observed i?-band magnitude reaches 
mj^'"'*=24. To determine the K-correction in the i?-band, 
we have redshifted the minimum reduced x^ spectrum 
model produced in our stellar mass calculations to mea- 
sure the i?-band magnitude at different distances. 



As the sensitivity of the X-ray telescope s hows energy 
and positional dependence (for details, see lYang et all 
l2004f ). to obtain a simple estimate of the solid angle at 
a give n hard X-ray flux, we have followed iBarger et al.l 
(J 2005l ) by comparing the observed number of X-ray ob- 
jects at different har d X-ray fl uxes with the av erage X-ray 
numb er counts from lCowie et al. (2003 ) and lYang et al.1 
()2004f ). The sample for comparing number counts is com- 
posed of all spectroscopically observed (not only those 
with secure redshifts) X-ray objects in the field of CDF- 
N, CDF-S and CLASXS and aU X-ray sources in the 
AEGIS field. Table [H lists the number of X-ray objects 
and the number of spectroscopic targets in all the fields. 
The spectroscopic targets are randomly selected for the 
CDF-S and CLASXS fields. Although the target selec- 
tion in the CDF-N field has a little bias, the majority of 
sources (439 out of 503) have been included in the spec- 
troscopic sample, making the effect of this bias on th e 
solid angle measurement neg ligible (iBarger et al.ll2005( l. 
Note that we have used the [Alexander et al.l ()2003[ ) X- 
ray catalog in the CDF-S due to its high X-ray positional 
accurac y so that our catalog contains ^30 less objects 
than in JBarger efaLl (|2005f) . The AEGIS field contains 
1318 objects and the redshift is obtained by matching the 
X-ray catalog to the DEEP2/AEGIS catalog whose tar- 
get selection depends on apparent magnitude and color. 
Therefore, for the objects in the AEGIS field, we have 
applied weights as shown below. Fig. [3] shows the ef- 
fective solid angle as a function of the hard X-ray flux. 
At the flux of 10"^'' erg s^^ cm^^, it is dominated by 
the CLASXS (0.4 square degree) and AEGIS fields (0.67 
square degree). 

We have followed the method in Willm er et al.l ()2006l ) 
to determine the galaxy weights, w, for the X-ray objects 
in the AEGIS field. These weights are used to account for 
the under-sampling of the photometric catalog (i.e. the 
fraction of objects with spectroscopic observations) and 
the redshift success rate of the spectroscopic catalog (i.e. 
the fraction of spectroscopic targets with secure redshift) 
of the DEEP2 observations. Briefly, we define a X-ray- 
optical photometric catalog as a ll AEGIS X.-ray o bjects 
that have optical counterparts in lCoil et al.l (|2004[ ). Here 
we assume that AEGIS X-ray objects without optical 
counterparts (about 25% of the whole X-ray sample) are 
not of interest for our study of low-mass systems at z < 
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Fig. 4. — The broad line AGN fraction in two rcdshift intervals 
of < ^ < 0.4 and 0.4 < z < 0.7. 

0.7, either at higher redshift or with steUar mass lower 
than our completeness limit (lO^'^ Mq). To demonstrate 
this, we found that 15% of the X-ray objects without op- 
tical counterparts have L^ > 10"*^ erg s~^ at z=0.7, i.e. 
satisfy the AGN definition. However, none of these ob- 
jects have absolute i?-band magnitude for mR = 24.75 at 
z=0.7 brighter than any of the observed low-mass AGN 
hosts with 9.7 < Log(M*/MQ) < 10.3, where m^ = 24.75 
is t he completeness cut of the DEEP2 photometry cata- 
log (|Coil et al.ll2"00l . For each AEGIS X-ray object with 
a secure spectroscopic redshift and probability of being 
a galaxy Pgai > 0.2, a data cube in the magnitude-color- 
color {R-{B-R)-{R-I)) space is defined. Then a proba- 
bility for being within the permitted redshift limit [0.1, 
1.4] is calculated for each object in the X-ray optical pho- 
tometric catalog within this data cube. The sum of all 
these probabilities within the data cube divided by the 
number of successful redshifts gives the weight uj, which 
is further corrected for the probability of being placed in 
the slit mask. 

As discussed in § [51 type 1 AGNs are excluded from the 
sample and thus a weight is applied to correct for their 
omission. The fraction of type 1 AGNs as a function 
of the X-ray luminosity is constructed in both redshift 
intervals as shown in Fig. [H which is almost the same as 
the result obtained bv Barger et al.l (|2005^ . The weight is 
calculated for each object based on the X-ray luminosity. 
As shown in Fig.[31 the corrections are small, as our AGN 
sample does not extend toward high X-ray luminosity 
where the broad line fraction is large. Therefore, our 
results throughout the paper do not change significantly 
if the broad-line AGN fraction does not apply. 

Finally, the mass function and X-ray luminosity func- 
tion are determined as: 



$(X)dX = ^L^/K„axdX, 



(3) 








22 24 26 
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where X is the stellar mass or X-ray luminosity, respec- 
tively, UJ is the galaxy weight, and Knax is the maximum 



Fig. 5. — The number distribution of all spectroscopic targets 
(open) and objects with secure spectroscopic redshift (hatched). 

volume for its detection. 

4.2. Incompleteness 

Our AGN sample where the 1/T^ax method applies are 
objects with secure spectroscopic redshifts, R <24, M* 
> 5 X 10^ Mq, and Ls^-kcV > 10^^ erg s'^. The criteria 
of secure spectroscopic redshifts does not introduce any 
significant incompleteness. Fig. [5] shows the distribution 
of R-band magnitude for the spectroscopically observed 
X-ray objects (open histogram) and the objects with se- 
cure spectroscopic redshifts (filled histogram). The red- 
shift success rate is 76% and 68% for spectroscopically 
observed objects at R <24 and 22 < R <24, respectively. 
Most of the spectro scopically failed objects have 22 < R 
<24. iBarger et al.l (2005) have shown that most objects 
with failed spectroscopic redshifts have photometric red- 
shifts larger tha n ^1. A similar result is found in DEEP2 
(jWillmer et al.l 2006). A rough estimate shows that all 
the spectroscopically failed objects contain roughly only 
two low-mass AGN hosts at z < 0.7, given that ~25% 
of the spectroscopically failed objects with 22 < R <24 
have a photometric redshift smaller than 0.7 and that 5% 
of the objects with secure redshift with 22 < R <24 are 
of low mass. 

The limiting i?-band magnitude of 24 is deep enough to 
sample most low-mass galaxies with M* > 5 x 10^ Mq 
out to a redshift of 0.7, the upperlimit where the spa- 
tial number density is measured. This is because 80% 
of AGNs with host Af* > 5 x 10® are detectable beyond 
a redshift of 0.7. The omission of a small fraction of 
low-mass hosts should not affect the comoving number 
density as they have been accounted for in the l/Knax 
method. To further demonstrate that the i? < 24 limit 
does not introduce any significant incompleteness even 
in the high redshift interval (0.4 < z < 0.7), we carried 
out a Monte-Carlo simulation to demonstrate that if the 
comoving number density is constant within 0.4 < z < 
0.7, the l/Vmax method can re-produce the intrinsic spa- 
tial density, regardless of the loss of a small fraction of 
faint R-band sources. Briefly, we randomly populated 
the volume with AGNs to a total number similar to that 
observed for 0.4 < z < 0.7. Their SEDs are assigned 
randomly to be one of those at 0.1 < z < 0.4 where the 
limiting i?-band magnitude of 24 should detect all low- 
mass galaxies with M* > 5x10® Mq. This is because 
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Fig. 6. — The stellar mass function of AGN hosts with i2-8fceV 
> 10''^ erg s^^ in the two redshift intervals of 0.1 < z < 0.4 and 0.4 
< z < 0.7. Dotted lines show the dividing mass between low-mass 
and massive hosts. The solid line in the upper panel shows the 
SDSS AGN result at z < 0.3 obtained by multiplying the fraction 
of galaxies hosting AGNs from Kauffmann et al. (2003a) with the 

f'axy mass function at 0.2 < z < 0.4 from Perez-Gonzalez et aT? 
08). The normalization of the AGN fraction of Kauffmann ct al. 
OSai ) has been decreased by a factor of 8 to match our normal- 
ization (a discussion about the difference betw een the SDSS and 
our low-redshift samples can be found in ^ 15.31 1. 



the iBruzual fc Chariot ()2003[ ) oldest simple stellar pop- 
ulations with solar metallicity and M^, = 5x 10^ Mq have 
mR=2A at z=0.4. Younger or lower metallicity galaxies 
will emit higher luminosities. We then applied the mea- 
sured z)i.™'* to the simulated galaxies and measured the 
1/Knax-based spatial number density. A thousand simu- 
lations show that there is no systematic difference in our 
estimate of the spatial density. For the incompleteness 
due to the X-ray detection threshold, since all AGNs (de- 
fined as i2-8kcV > 10''^ "^rg s~^) in low-mass hosts with 
M* > 5 X 10^ have 4™* > 0.7, the 1/V;„ax method should 
recover the real spatial number density. 

5. RESULTS 
5.1. Comoving Number Density and AGN Fraction 
Fig. [6] shows the spatial number density of AGNs with 
-^2-8koV ^ l'^^^ ^^S ^~^ ^s ^ function of AGN host mass 
in the two redshift intervals of 0.1 < z < 0.4 and 0.4 < z 
< 0.7. Dotted lines show the dividing mass between low- 
mass and massive hosts. In both redshift intervals, the 
AGN number density depends on the host mass, peak- 
ing at an intermediate mass range and decreasing toward 
higher and lower mass. Furthermore, the AGN host mass 
function peaks at a higher mass in the higher redshift 
interval. Our low redshift interval suffers from low num- 
ber statistics with on average ^7 objects per mass bin, 
compared to ^20 objects per mass bin in high redshift 
interval. However, the peak around a stellar mass of 
-[^Qio.3 _ ]^q10-8 ]y[(^ is quite possibly real. For example, 
iGreene fc Hd (|2007b( l have measured the local BH mass 
function of type 1 SDSS AGNs at z < 0.3. Although 
suffering from an incompleteness problem, there appears 
to be a turnover of their BH mass function around a BH 



mass of 10^ Mq, which corr esponds on avera ge to a host 
stellar mass of lO^^^ Mq (|Best et al.l 120051 ). The sta- 
tistical significance for the low-redshift AGN host mass 
function peaking in the second mass bin was estimated 
using a Monte-Carlo simulation. Basically, we perturbed 
the mass function at each mass bin ten thousand times 
assuming a normal distribution with the 1-cr deviation 
equal to the measured error. The probabilities for the 
peak at first, second and third bins are 10%, 80% and 
10%, respectively, which indicates that the low-redshift 
AGN host mass function most likely peaks in the second 
mass bin. 

To further demonstrate that the AGN host mass func- 
tion peaks at a higher mass in our high redshift inter- 
val, the solid line in the upper panel of Fig. [6] shows 
the SDSS AGN result at z < 0.3 obtained by mul- 
tiplying the frac t ion of galaxies hosting AGNs from 
iKauffmann et all (l2003al) with the galaxy m ass function 
at 0.2 < z < 0.4 from lPerez-Gonzalez etal. (200 g). The 
galaxy mass function obtained by iPerez-Gonzalez et al.l 
Hops) is based on a rest-frame near-IR selected galaxy 
sample and is generally consistent with results of other 
studies (any differences are with in 0.3 dex (see their 
Fig. 4) ). At z < 0.7, the result of iPerez-Gonzalez et al.l 
(|2008?) is complete down to stellar mass of 10^ M0, deep 
enough for our purpose. The norma lization of the AGN 
fraction of IKauffmann et all (|2003al ) has been decreased 
by a factor of 8 to match our normalization (a discus- 
sion about the difference between the SDSS and our low- 
redshift samples can be found in § 15. 3|) . With much 
higher number statistics, the SDSS result shows the AGN 
host mass function peaks around 10^°'^ M©. We then 
carried out a simulation by assuming the 0.4 < z < 0.7 
AGN mass function actually follows the SDSS result at z 

< 0.3. For each object at 0.4 < z < 0.7, a host mass is as- 
signed randomly with a relative probability that follows 
the SDSS result. The range of simulated stellar mass is 
from 10^ to 10^^-^ Mq. We also assumed that the total 
probability in this mass range is equal to 1. Combining 
the simulated stellar mass and the observed zi.™'* and 

XI ay 

^hmit nieasured in § |4l we can construct the AGN mass 
function at 0.4 < z < 0.7 using the same mass bins as 
for the observed data. Ten thousand simulations indi- 
cate that the probability is 99.6% for the mass function 
to peak in one of the two lowest mass bins. That is, the 
simulated host mass function at 0.4 < z < 0.7 (using 
the SDSS low-z data as the basis of the simulated data) 
would hkely (99.6%) peak in the two low-mass bins (i.e. 

< 10^^ M0), while the observations at 0.4 < z < 0.7 find a 
significant peak at higher mass (> lO'^^ M0). This result 
indicates that the tendency for the AGN mass function 
in the high redshift interval to peak at higher mass is 
significant (99.6%). 

The simulations shown above test the pure number 
statistics. In order to account for the possible selection 
bias that the low-mass galaxies harbor faint AGNs due 
to the mild correlation between the galaxy stellar mass 
and the central BH mass coupled with a given Eddington 
ratio distribution, we first created a set of stellar masses 
with the relative probability following the galaxy stellar 
mass function from Perez- Gonzalez et al. (2008). These 
stellar masses were then converted to the BH mass dis- 
tribution using the average ratios of BH masses to galaxy 
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Fig. 7. — The fraction of galaxies hosting active nuclei as a 
function of host stellar mass in the two redshift intervals of 0.1 < 
z < 0.4 and 0.4 < z < 0.7. Dotted lines show the dividing mass 
between low-mass and massive hosts. The solid line in the upper 
panel shows the trend of the SDSS powerful AGN fraction with the 
host stellar mass from Kauffmann et al. (2003a) whose normaliza- 
tion has been decreased by a factor of 8 (a discussion about the 
difference between the SDSS and our low-rcdshift samples can be 
found in ij 15.311 . 



stellar masses as a function o f stell ar masses and associ- 
ated errors from iBest et al.l (|2005[ ). For the Eddington 
ratio between [10~^, 1], we assumed a probability distri- 
bution of P(Lboi/iEdd) d(Lbol/iEdd) OC l/(Lbol/iEdd)", 
where n is a free parameter. By assuming Lboi/^x = 20, 
the simulated AGN sample is then defined as objects with 
Lx > 10^^ erg s~^ and the host mass function of this sam- 
ple can be constructed. We simulated about ten times for 
each n value which is given in a set of discrete numbers 
starting at zero and increasing in a step of 0.5. At n—2, 
the simulated AGN mass function shows a similar trend 
to our measured mass function in the redshift interval of 
0.4 < z < 0.7. The peak number density is about four 
times larger than the value around the lowest measured 
mass (10^" Mq). This result shows that when the Ed- 
dington ratio distribution is strongly biased toward a low 
value, the observed turnover of the AGN mass function 
is due to the selection bias. However, the corresponding 
simulated fractions of galaxies hosting AGNs are an order 
of magnitude lower than the observed ones. This implies 
the turnover of the AGN mass function most likely in- 
dicates a real decrease of the number density of galaxies 
with active BHs. 

Compared to the galaxy mass function, the unique 
feature of these two AGN mass functions is that they 
do not increase monotonically with decreasing mass. A 
physical parameter more related to the AGN cycle and 
BHOF is the fraction of galaxies hosting active nuclei. 



which is defined by dividing the AGN host mass func- 
tion by the galaxy mass function at a given stellar mass, 
as shown in Fig. El The 0.2 < z < 0.4 and 0.4 < z 
< 0.6 galaxy mass functions from iPerez-Gonzalez et al.l 
()2008f ) are used for our low and high redshift intervals, 
respectively. Fig. [7] shows that more massive galaxies 
host active nuclei more frequently in both redshift in- 
tervals. In the local universe, it has been found that 
the AGN fraction depends on the host stellar mass 



jKauffmann et al."2003a UGallo et "aI1l2007HDecia,rH et al 



2007; Sivakoff et al. 2001. ho st Hubble t vpe (iHo et al. 



1297) and BH mass (jHeckman et al.ll2004l : iGreene fc He 
2007b) . Note that these three parameters are roughly 



correlated with each other. In general, the fraction har- 
boring active nuclei is lower for lower-mass local systems. 
Our result indicates that such a dependence also exists 
at high redshift. However, there are dif ferences in the be- 
havior between low- and high-redshift. I Kauffmann et all 
()2003a[ ) show their z < 0.3 SDSS AGN fraction increases 
with stellar mass from low-mass galaxies up to galaxies 
with 10^^ Mq and then drops quickly toward higher mass 
galaxies, which is shown as the solid line in the upper 
panel of Fig.[7l Our result for 0.1 < 2: < 0.4 is consistent 
with their work, particularly since it indicates a leveling 
out of the AGN fraction near 10^^ M0. However, our 
sample contains too few hosts above this limit to con- 
firm the drop toward even higher masses. At 0.4<z<0.7, 
our study shows that the AGN fraction increases all the 
way from low-mass galaxies to massive galaxies with M* 
— 10^^ Mq. Note that our AGN fraction for massive 
galaxies (M, > 10^^ Mq) is consiste nt with other stud- 
ies of X-ray AGNs at this redshift (jBundv et all l2007t 
lAlonso-Herrero et al.ll2008[ ). while showing that this lack 
of turn-over continues up to higher mass. The difference 
between our high redshift interval and the SDSS local 
universe may be explaine d by the cosmic "d ownsizing" 
evolution of MBHs (e.g. iBarger et al]l2005( l. Massive 
BHs in massive galaxies have accumulated their masses 
at high redshift and become quiescent in the local uni- 
verse. 
Note that our AGNs are defined as those with i2-8koV 

> lO"*^ erg s^^. It is obvious that the AGN number 
density or AGN fraction depends on the depth of sur- 
veys, or, equivalently, the limiting Eddington accretion 
ratio for a given BH mass. As shown in § |3l the limit- 
ing Eddington ratio is 0.02 for our lowest stellar mass 
of 10^'^ Mq. For a galaxy with stellar mass of 10^^ 
Mq, the limiting Eddington ratio is 0.001. To have a 
sense of the missed fraction of AGNs with lower Ed- 
dington ratios, we use d the cur rently most complete lo- 
cal AGN sample from lHo et al.l (1995) as a comparison. 
Our sample of AGNs should miss a significant fraction 
of high-ionization Seyfert galaxies, as local Seyferts have 
a median Eddington ratio of 1.3x10"'' (|Ho| [2008). We 
may miss all low-ionization nuclear emission-line regions 
(LINERs), as none of the local LINERs has an Edding- 
ton ratio > O.OOl (Ho 2008). OveraU, only 5% and 10% 
of all local AGNs have Eddington ratios larger than 0.01 
and 0.001, respectively (Ho 2008|). Therefore, the frac- 
tion of galaxies with active BHs at high redshift may be 
much larger than our result, which is limited to L^2-8kc'v 

> W^^ erg s~^. Identifying these additional AGNs is dif- 
ficult because of the overlap of their X-ray luminosities 
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Fig. 8. — The rest-frame 2-8 KeV X-ray luminosity function 
of AGNs in low-mass hosts, and of all AGNs in the two redshift 
intervals of O.K z <0. 4 and 0.4< z <0.7. T he two solid lines show 
the XLF obtained by IBarger jFan l(200l ). based on the CDFN, 
CDFS and CLASXS fields. 



with those resulting from star formation (iRanalh et al.l 
[2003). However, it is still important to understand the 
fraction of galaxies hosting AGNs with relatively high 
Eddington ratios and their trends with host masses, as 
these AGNs may be in the phase when feedback from 
MBHs is most important and also the main stage of BH 
growth through accretion. 

5.2. X-ray Luminosity Function of AGNs in Low-Mass 

Hosts 

Fig. [8] shows the X-ray luminosity functions of AGNs 
in low-mass hosts (diamonds) and all AGNs (asterisks) 
in the two redshift intervals of 0.1 < z <0.4 and 0.4< 
z <0.7. Tw o solid lines show the XLFs obtained by 
IBarger et al] ([2005h . based on the CDFN, CDFS and 
CLASXS fields. The XLFs of our AGN host sample with 
i? <24 closely match their results for the higher redshift 
interval. In the low redshift interval, there is about a fac- 
tor of two difference in the highest luminosity bin. This 
is most likely caused by the cosmic variance due to the 
small comoving volume at low redshift. 

The total X-ray energy density of AGNs in low- 
mass hosts can be measured by integrating the XLF 
where there are data points. The measurements give 
(2.2±1.4)xl03^ and {lA±OA)xlO^^ erg s'^ Mpc^^ at 
0.1 < z < 0.4 and 0.4 < z < 0.7, respectively. Similarly, 
the X-ray energy density of all AGNs can be obtained. 
The contribution of AGNs in low-mass hosts to the to- 
tal X-ray energy density is (11±8)% and (14±5)% in the 
lower and higher redshift intervals, respectively. There- 
fore, these AGNs are energetically important to the cos- 
mic X-ray background. 
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Fig. 9. — The redshift evolution of AGNs with host stellar mass 
9.7 < logM*/MQ < 10.3. From top to bottom: the comoving 
number density; the fraction of low-mass galaxies hosting active 
nuclei; and the comoving X-ray energy density. 

5.3. Redshift Evolution of AGNs in Low-Mass Hosts 

As shown in Fig. [51 AGNs in low-mass hosts show 
strong redshift evolution. The number density increases 
with redshift following $dagn = 10-5°±°-4 (l-Hz)30±2-5 
Mpc~'^ log(M0)~^, the fraction of low-mass galaxies with 
AGNs evolves as /agn = io-2.9±o.4(;l_^_^)3.6±2.5 ^^^ ^^le 
energy density follows pdagn = io36.5±o.6(^_,_2)8.5±3.2 
erg s~^ Mpc~'^ in the redshift range of 0<z<0.7. Note 
that these three quantities are related to each other to 
some degree. 

To evaluate our result for cosmic evolution of AGNs 
in low-mass galaxies, we compared it at low-redshift 
to those based on SDSS optical-emission-line-selected 
AGNs. Greene & Ho (2007b) show that the spatial num- 
ber density of BHs between 10^- ^ and lO'''-^ Mp hosted 
by SDSS narrow-hne AGNs from lHeckman et all (|2004D 
is around 6x10"^ Mpc"^ logM^^'j (see their Fig. 10). As- 
suming the host galaxies of these BHs have M^ ^ 10^° 
Mq, this number density is three times higher than the 
prediction of -(2.0±1.5)xl0"^ Mpc^^ logM"! at a red- 
shift of 0.15 based on the number density evolution of 
our AGNs in low-mass hosts. Note that our result has 
been corrected for broad-line AGN. However, the cor- 
rection is quite small ('^10%) as shown in § 14.11 and 
Fig. [H The difference could be partly caused by our 
low number s tatistics. A lthough the limiting Eddington 
ratios of the I Greene fc Ho (2007b)'s AGNs are similar 
to ours (^ 0.01), the different selection algorithms (op- 
tical emission diagnostics vs. X-ray emission) may be 
another reason for the discrepancy. The optical emis- 



sion diagnostics miss AGNs when the contrast of narrow 
Hne emission to host galaxy hght is low. At low red- 
shift, the slit for the spectroscopic observations covers 
only part of the host galaxies so that host galaxy light 
dilution is not quite so severe. iKauffmann et all (J2Q033) 
show that the SDSS optical emission line diagnostics can 
recover the AGNs with Eddington ratios below 0.01 for 
host galaxies with M, ~ 10"'^° Mq, if the emission lines 
are not strongly extincted. The X-ray emission method 
misses AGNs whose X-ray emission is obscured so that 
their observed X-ray luminosity falls below our thresh- 
old for AGN definition. It is obvious that Compton- 
thick AGNs {Nh > 10^"* cm"^) ^jn ^g missed in the 
2-10 keV X-ray survey. However, if AGNs in low-mass 
hosts have intrinsically weaker X-ray emission compared 
to more massive host AGNs, the heavily-extincted but 
Compton-thin AGNs (lO^^ ^^-2 < at^ < io24 cm'^) 
may still be missed. Our X-ray spectral fits of AGNs in 
low-mass hosts with X-ray counts > 100 (Shi et al. 2008, 
in preparation) have shown tha t only ~10% o f these 
AGNs have Nm > lO^^ ,,^^-2 [Heckman et"all (|200l 
did show that the emission-line-selected AGN sample is 
more complete than the X-ray-selected one at low red- 
shift and that the SDSS AGN spatial density is about 
three times higher than X-ray-selected AGNs. Although 
at low redshift optical emission line diagnostics are more 
complete, X-ray selection identifies more AGNs at high 
redshift. X-ray selection also provides more uniform se- 
lection criteria (e.g., less affected by dilution by the host 
galaxy extended emission). As a summary, the prediction 
at low-redshift based on our result of cosmic evolution of 
AGNs in low-mass galaxies is generally consistent with 
the SDSS result, with an offset caused by different selec- 
tion biases between X-ray-emission-selected AGNs and 
optical-emission-line-selected ones. 

6. DISCUSSION: BLACK HOLE OCCUPATION FRACTION 
IN LOW-MASS GALAXIES 

6.1. Fraction of Galaxies Hosting AGNs 

§ 15.11 shows that the AGN mass function peaks at an 
intermediate mass and decreases toward the low-mass 
regime, in contrast to the monotonic increase of the 
galaxy mass function. § 15.11 also indicates that the frac- 
tion of galaxies hosting active nuclei decreases with de- 
creasing host stellar mass. Theoretically, the fraction of 
galaxies hosting AGNs can be estimated by 

/agn = /bhof7^agn (4) 

where /bhof is the black hole occupation function, i.e., 
the fraction of galaxies hosting MBHs at their centers, 
7 is the fractional AGN trigger rate, i.e., the fraction 
of MBHs becomes active per unit time, and ^agn is the 
duration of a nuclear activity episode with Lx>lO^^ erg 
s^^. The trend of AGN fractions with host mass can be 
caused by the host mass dependence of any one of the 
three factors, /bhof, 7 and ^agn- Different models for 
MBH seed formation in the early universe predict differ - 
ent mass dependences of /bhof (jVolonteri et al.l[2008f ). 
^AGN is most likely mass dependent, as the lower mass 
systems require larger Eddington ratios to be brighter 
than Lx=W^^ erg s~^. 

To get some idea of /bhof of low-mass galaxies, we 
make a simple assumption that each major merger trig- 
gers one-time nuclear activity and thus 7 = the merging 



rate. Although the measurement of the merging rate is 
subject to various uncertainties (e.g. lCassata et al1l2005t 
IShi et al.ll2006tlLotz et aI]|2008D . ft should be around 0.1- 
0.2 Gyr~^ for massive galaxies. Simply assuming 7 = 0.1 
Gyr~^ for low-mass galaxies, we have : 



O.OSGyrO.l ,^ ^ ^ ^, 

/bhof = — — — ; (0.1 < z< 0.4) 



/ 



BHOF — 



^AGN 7 

O.OGGyrO.l 



; (0.4 <z< 0.7) 



(5) 



(6) 



^AGN 7 

In this simplified case, as long as the duration of one 
episodic nuclear active phase is not long (<0.06 Gyr), all 
low-mass galaxies with 9.8 < logM* < 10.3 at z < ~1 
should harbor BHs at their centers. 

6.2. Accreted BH Mass by AGNs in Low-Mass Galaxies 

The amount of BH mass accreted by AGNs in low- 
mass galaxies since z — I can be measured by inte- 
grating the cosmic evolution of the X-ray energy den- 
sity of these AGNs obtained in § 15.31 This mass 
should provide a stringent lower limit to the BHOF 
in local low-mass galaxies. An average X-ray to bolo- 
metric luminosity correction of 19 is obtained by us- 
ing Lboi/L2-inkeV = 17(L2-iokev/10'*^ erg 5-^)°-*^ 
(IShankar et al.|[20M ). L2-8kev/i2-iokoV=0.86 and the 
mean X-ray luminosity of 7x10''^ erg s~^. Assuming 
the mass-to-radiation conversion efficiency e=0.1, the 
total accreted black hole mass in galaxies with 9.7 < 
log(M*/Mo) < 10.3 since z = 1 is 3.9(±0.9)xl03 Mq 
Mpc~3. 

This accreted BH mass must be hosted in local low- 
mass galaxies. The corresponding total galaxy mass can 
be estimated using the bulge-BH relation and bulge-to- 
disk ratio. The lower limit to the local BHOF is then the 
ratio of the total mass of galaxies hosting the accreted 
BHs to the total local mass in low-mass galaxies. In 
practice, we can alternatively assume all local low-mass 
galaxies host BHs. The lower limit of the BHOF is then 
the ratio of the accreted BH mass to the assumed total 
BH mass hosted by all local low-mass galaxies. 

The BH mass function can be determined from the 
galaxy velocity dispersion function using the relationship 
between the velocity dispersion and BH mass (Mbh-o')- 
However, it is relatively difficult to measure the veloc- 
ity dispersion, resulting in incompleteness of the BH 
mass function at the low mass end. Alternatively, it 
can be derived from the bulge luminosity function us- 
ing the relationship between the bulge luminosity and 
the BH mass (LbuIsc-Mbh)- In this case, the bulge-to- 
disk ratios of galaxies with different morphologies are 
required to derive the bulge luminosity function from 
the galaxy total luminosity function. In practice, the 
two methods are actually employed at the same time to 
complement each other. Different studies produce a rel- 
atively consistent result for the total local BH mass den- 
sity, which is around 4.5 x 10^ Mq Mpc~'^ with uncer- 
tainty < 2.0x10^ M.T^ Mpc"^ (lAiler fcRichstoiIg 120021: 
IShankar et al.ll2004 [Marconi et al.ll2004f ). Therefore, the 
accreted BH mass by AGNs in low-mass hosts since z — 
1 is only a small fraction (0.8%) of the total local BH 
mass density. 

To estimate the BHOF in local low-mass galaxies 
using the accreted BH mass during their AGN phase 
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since z ~ 1, we first need to determine the local 
successor of our low-mass AGN hosts. Given the cosmic 
evolution of the stellar mass densit y of galaxies with dif- 
ferent stellar masses provided by iPerez-Gonzalez et alj 
(|2008[ ). the 9.7 < log(M,/MQ) < 10.3 interval at 
median redshift oi z = 0.5 corresponds roughly to 
9.85< log(MH,/M0) < 10.45 at z = 0. To estimate 
the total BH mass associated with local galaxies with 
9.85 < log(Af*/Af0) < 10.45, we measured the BH 
mass hosted by the early-type and late-type galaxies 
separately, as they have different bulge to total light 
ratios, f^'^^J^ = 0.8 5±0.05 and /^=^'i°^ = 0.30±0.05 
(jShankar et all \2004) . Assuming that log{MBH/MQ) 

= (8^0±gjI))-H(1.12±0.06)log(Afbuigo/10"Mp) 

(jHaring fc Rixl l2004f ). early and late type galaxies 
with 9.85 < log(MH./MQ) < 10.45 harbor BHs with 6.85 

< log{MBH/MQ) < 7.50 and 6.15 < log fMBg/M^) 

< 6.8 0, respectively. By fitting Fig. 5 of IShankar et al.l 
(|2004t) with the formula 



BH 



Mbh s 

^M )' 

Mbh 



,/3l 



(7) 



we obtained <i>*gjj = 5.2(±0.2)xl0"3 Mpc'^, Af^^ = 
8.8(±2.0)xlOS Mq, a = -0.47(±0.03) and f3 = 
0.39(±0.01) for early-type galaxies, and ^*bh = 
1.9(±0.2)xl0~2 Mpc-3, M^iH = 3.3(±1.7)xl0'5 M©, a 
= -0.57(±0.09) and P = 0.34(±0.02) for late-type galax- 
ies. Integrating the above equation over the range of in- 
terest, the local early- and late-type galaxies with 9.85< 
log(M,/M0) < 10.45 harbor total BH masses of 2.1 x lO'' 
and 1.2x10^ Mq Mpc~'^, respectively. Therefore, the ac- 
creted BH mass by AGNs in low-mass hosts from z = 1 
contributes 12% of the local BH mass hosted in galaxies 
with 9.85< log(M,/MQ) < 10.45. 

Three main uncertainties associated with this percent- 
age include the accreted BH mass, the stellar mass of 
local counterparts of our high redshift low-mass hosts 
and the local BH mass function. As shown above the 
measurement error of the accreted BH mass is about 
20%. A larger uncertainty for this accreted BH mass 
is from the uncertain estimate of AGN accretion in low- 
mass galaxies at 0.7 < z < 1. The current assumption 
is that it follows the trend at z < 0.7. If the energy 
density of AGNs in the low-mass host is flat at z > 0.7, 
the accreted BH mass will decrease by about a factor 
of 2. To estimate the second uncertainty, we can as- 
sume an extreme case where all our AGNs in low-mass 
hosts are at z = 1. In this case, the local counterparts 
have 10. 1< log(M*/M0) < 10.7 and the accreted BH 
mass contribution is 6%. The uncertainty in the local 
BH mass functio n is about a factor of 3. from comparing 
different studies (lAller fc Richstonell2002l : [Shankar et all 
120041: iMarconi et al.ll2004D . Therefore, the final uncer- 
tainty could be a factor of 4, estimated by adding the 
above errors quadratically. 

A more important factor, resulting in severe underes- 
timate of the accreted BH mass by AGNs, is the X-ray 
obscuration correction for the observed X-ray flux plus 
the significant number of missed heavily-obscured AGNs 
in current X-ray surveys. This problem is even worse for 
AGNs in low-mass hosts, as (1) they on average have low 
intrinsic X-ray luminosities, with many falling below our 
selection criteria of Lx>W^^ erg s^^ and (2) the struc- 



ture of the circumnuclear material may vary with the 
luminosity in a way that causes lower luminosity objects 
to be more obscured, as implied by the decrease in the 
fraction of type 2 o bjects with increasing luminosity (e.g. 
'Barge r et al.ll2005l ) . We calculated the degree of under- 
estimation in two ways. First, we measured the accreted 
BH mass of 3.4x10^^ Mq Mpc"^ by all X-ray-detected 
AGNs from z = 1 by integrating the X-ray energ y den- 
sity e volution of all AGNs (see § 15. 2p . Shankar et al.l 
([2004") obtained total BH mass of -1.3x10^ M0 Mpc'^ 
accreted by all AGNs from z = 1 with a correction for 
obscuration and accounting for missing Compton-thick 
objects. Therefore, the accreted BH mass by AGNs in 
low-mass galaxies is underestimated by a factor of 4, by 
assuming that the obscuration of these AGNs is simi- 
lar to that of all AGNs. The second way is to compare 
the distribution of the HI column density of observed 
AGNs to the intrinsic distribution to estimate the frac- 
tion of missing AGNs. Our X-ray spectral fits of AGNs 
in low-mass galaxies with X-ray counts > 100 (Shi et al. 
2008, in preparation) have shown that only 10% of such 
AGNs have HI column density A^m > 10^^ cm"^. The 
study of local Seyfert galaxies sho ws that 75% of Sey fert 
2 galaxies have Nm > lO^^ cm'^ (jRisaliti et al.lll999h . If 
AGNs in low-mass galaxies have structures of circumnu- 
clear material similar to those in Seyfert galaxies, 75% of 
the low-mass galaxies should be missed and the accreted 
BH mass by all such AGNs is underestimated by a factor 
of 4, similar to the value estimated by the first method. 

Therefore, the accreted BH mass by all AGNs in low- 
mass hosts including missed heavily-obscured ones from 
z = 1, probably contribute ^ 50% of the local BH mass 
hosted in galaxies with 9.85< log(M*/M0) < 10.45. As 
discussed at the beginning of this section, this percent- 
age is obtained by assuming all such galaxies harbor BHs 
at their centers. A percentage of 50% can also be inter- 
preted as the lower limit of the BHOF in local galaxies 
with 9.85 < log(M*/M0) < 10.45, i.e., at least 1 out of 
2 local galaxies with 9.85 < log(M*/M0) < 10.45 har- 
bor BHs. Accretion at higher redshift (z > 1), at low 
Eddington ratios, or the mass growth through the BH 
merging processes may provide the remaining BH mass. 

The above results strongly favor a scenario where the 
nuclear activity is a necessary ingredient during low- 
mass galaxy evolution. The fraction of low-mass galaxies 
hosting active nuclei (see § 16. ip and the amount of the 
accreted black hole mass in low-mass host AGNs (see 
this section) implies that a significant fraction of low- 
mass galaxies in the local universe harbor black holes at 
their centers and these black hole masses are assembled 
through the active accretion phase {Lx > 10^^ erg s""'^) 
at z < 1.0. 

7. CONCLUSIONS 

We have studied a sample of X-ray AGNs in low-mass 
host galaxies with stellar mass of 5x10^ M0 < M^, < 
2x10^" M0 out to z '^1. Our main conclusions are: 

(1) By including AGNs in more massive host galaxies, 
we have constructed the stellar mass function of AGN 
host galaxies extending down to the low-mass regime in 
two redshift intervals of 0.1 < z < 0.4 and 0.4 < z < 0.7. 
We have found the AGN host stellar mass function peaks 
at an intermediate mass range, with the peak shifting 
toward higher mass at higher redshift. 



11 



(2) By comparing to AGNs in more massive hosts, we 
have found that the fraction of galaxies hosting active 
nuclei depends strongly on the host mass, increasing from 
low-mass galaxies continuously to the most massive ones. 
The fraction of low-mass galaxies hosting active nuclei 
suggests that a large fraction of such galaxies at < z 
< 1 harbor black holes at their centers as long as the 
low-mass host AGN lifetime with Lx > W^^ erg s~^ is 
not long (< 0.06 Gyr). 

(3) The X-ray luminosity functions of AGNs in low- 
mass hosts have been constructed in two redshift inter- 
vals of 0.1 < z < 0.4 and 0.4 < z < 0.7. AGNs in low-mass 
hosts contribute ~10% of the X-ray energy density of all 
the AGNs in < z < 1, indicating that such AGNs make 
an energetically significant contribution to the cosmic X- 
ray background. 

(4) AGNs in low-mass hosts show strong redshift evo- 
lution in their comoving number density, the fraction of 



such galaxies with active nuclei and the comoving radia- 
tion energy density. By integrating the X-ray luminosity 
function of these AGNs over the redshift range of < z 
< 1, the accreted black hole mass in galaxies with 5x10^ 
Mq < M, < 2x101° Mq is (3.9±0.9)x 10^ Mq Mpc'^. 
This number gives a strong lower limit of 12% to the 
fraction of local low-mass galaxies harboring black holes, 
which may be much higher ( > 50%) if the dusty torus 
of the AGNs in low-mass galaxies has similar structure 
to or is more opaque than that of AGNs in massive host 
galaxies. 



We thank the anonymous referee for detailed com- 
ments. Support for this work was provided by NASA 
through contract 1255094 issued by JPL/ California In- 
stitute of Technology. 
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TABLE 1 
X-RAY Fields 
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FIELD 



(1) 



Area 

dcg^ 

(2) 



TTtlimit 

erg s ^ cm 



Rcf Optical photometry 



(3) 



(4) 



(5) 



Rcf 



(6) 



Secure redshift Ref 

(7) (8) 



AEGIS 


0.67 


S.OxlO-if' 


1 


B, R, I, Ks 


2,3 


CDF-N 


0.12 


1.4x10-1'^ 


5 


U, B, V, R, I, z, HK 


6 


CDF-S 


0.11 


2.8x10-18 


8 


U, B, V, R, I, 914nm 


9 


CLASXS 


0.4 


3.0x10-15 


11 


B, V, R, I, z 


12 


XMMLSS 


1 


5.8x10-15 


13 


B, V, R, I, J, K 


14,15,16,17 



zquality=3, 4 4 

all z without label 's' 6, 7 

qual=2 10 

all listed z are secure 12 

q_z=3, 4, 13, 14, 23, 24 18 



Note. — Col.(l): The field name. Col. (2): The area of the field size in square degree. Col. (3): The limiting X-ray 
flux in 2-8 keV. Col. (4): The reference for the X-ray data. Col. (5): The available optical/near-IR photometry. Col. (6): 
The reference for the optical/near-IR photometry. Col. (7): The definition of the secure spectroscopic redshift in each 

field. Col. (8): T he reference for the s pectroscopic redsh ift. 

References - (l)|Nandra et al.l |[2005h : JLaird et al. | II 2008I): (2)ICoil et al.l (20041): f SllBundy et al.l ll2006l): (4'llDavis et al.| 
( 20031 . [2007tl: (5 ) Alexanc to^et al.l I I2003D: (e) Barger et al. (2003); (7) Barg er et al.l | |200S ); ( 8) Alexander et al.l 
jgOOsh : (9) IWolf et al. (200l) : ClO^ ISzokoly et al ( 2004): (11) Yang et al. I2004h : (1 2'1 IStefTen et al. ( 2004) ; (13) 
Chiappetti et al. (2005): (14) Le Fcvrc et all J2004|) ; (15) IMcCracken et al.l | |2003|) : (16') IRadovich et al.l (120041) : (17) 



lovino et al., (,20p5i); (18) ,Le Fcvre et al., (.20051) 



TABLE 2 

The number of X-ray objects in all fields 



FIELD CDF-N 


CDF-S 


CLASXS 


AEGIS 


XMMLSS 


Total 503 
Spec-Observed 439 


326 
210 


525 

422 


1318 
357 


286 
23 



TABLE 3 
The Parameters of Stellar Synthesis Models. 



parameters 



value 



Simple stellar populations 

Metallicity 

visual extinction Tu 

e-folding time r for exponential star-formation history 

fraction of ejected gas to be recycled e 

galaxy age 



Chabricr (2003) IMF and Padova 1994 evolutionary tracks 

0.0001, 0.0004, 0.004, 0.008, 0.02 (Zq), 0.05 

[0.0, 4.0] with a step of 0.5 in logarithm 

[0.05, 8.91] Gyr with a step of 0.25 in logarithm, 100 Gyr 

0.001, 0.01, 0.1, 1 

[0.001, 15.85] Gyr with a step of 0.1 in logarithm 



TABLE 4 
Sample of AGNs in Low-Mass Host Galaxies 



(1) 



(2) 



RAx, DECx 

(3) 



D. 



opt 



(4) 



Phtomctry 

[in AB system] 

(5) 



[10- 



/x_ 
^crgs -'cm ^1 



(6) 



AEGIS123 


0.92 


14 15 39.2 +52 08 49.6 


0.98 


AEGIS232 


0.81 


14 16 16.8 +52 21 35.6 


1.37 


AEGIS375 


0.48 


14 17 24.6 +52 30 25.0 


0.31 


AEGIS555 


0.84 


14 18 25.2 +52 49 20.8 


0.77 


AEGIS655 


0.71 


14 19 06.5 +52 38 55.8 


1.29 


AEGIS723 


0.66 


14 19 30.8 +52 56 17.3 


0.42 


AEGIS795 


0.46 


14 20 01.4 +52 53 10.7 


0.96 


AEGIS969 


0.91 


14 20 59.8 +52 56 04.3 


0.92 


AEGIS1008 


0.43 


14 21 15.9 +53 19 48.6 


0.38 


AEGIS1159 


1.24 


14 22 14.8 +53 23 54.3 


1.90 


AEGIS1303 


i.on 


14 23 26.1 +53 30 03.7 


0.79 


0DF-N5 


n.56 


12 35 21.3 +62 16 28.1 


0.05 


CDF-N83 


n.46 


12 36 08.2 +62 15 53.1 


1.52 


CDF-N191 


n.56 


12 36 35.9 +62 07 07.7 


1.43 


CDF-N194 


n.56 


12 36 36.7 +62 11 56.0 


0.34 


CDF-N267 


n.4n 


12 36 51.7 +62 12 21.4 


0.97 


CDF-N301 


n.29 


12 36 58.7 +62 04 02.4 


0.62 


CDF-N441 


n.63 


12 37 36.0 +62 18 05.9 


0.00 


CDF-S241 


n.68 


03 32 39.1 -27 44 39.1 


1.62 


OLASXS42 


0.49 


10 31 54.9 +57 45 20.9 


0.69 


CLASXS131 


0.39 


10 32 42.6 +57 56 20.8 


0.72 


CLASXS205 


0.68 


10 33 18.1 +57 26 01.4 


0.65 


CLASXS231 


1.38 


10 33 29.2 +57 47 08.1 


0.56 


CLASXS243 


0.32 


10 33 34.1 +57 56 01.9 


0.18 


CLASXS286 


0.29 


10 33 53.2 +57 32 41.0 


0.48 


CLASXS322 


0.20 


10 34 06.6 +57 56 07.3 


0.26 


CLASXS329 


0.37 


10 34 09.5 +57 29 53.8 


0.30 


CLASXS373 


0.33 


10 34 29.7 +57 50 58.2 


0.69 


CLASXS441 


0.39 


10 34 56.2 +57 47 24.5 


0.85 


CLASXS448 


0.62 


10 34 57.9 +57 37 56.1 


0.42 


CLASXS517 


0.62 


10 35 51.0 +57 43 33.0 


0.28 


CLASXS522 


0.51 


10 36 04.2 +57 47 48.3 


0.17 



24.35 


23.95 


23.71 


24.00 


23.31 


22.60 


20.23 


19.99 


19.66 


24.48 


24.06 


23.38 


23.52 


22.65 


22.18 


25.99 


23.74 


22.76 


23.85 


21.91 


21.32 


23.06 


22.71 


22.00 


22.27 


21.23 


20.81 


24.02 


23.94 


23.35 


22.54 


21.89 


21.64 


23.20 


22.82 


22.42 


23.50 


23.22 


22.62 


25.40 


25.52 


25.22 


24.30 


23.42 


23.12 


24.20 


23.52 


22.62 


22.20 


21.52 


21.62 


25.20 


25.02 


24.52 


25.84 


25.32 


25.44 


24.90 


24.70 


24.10 


23.30 


23.90 


23.50 


21.50 


22.00 


21.70 


23.40 


23.60 


23.20 


21.90 


21.20 


20.50 


21.00 


20.90 


20.20 


21.50 


21.00 


20.50 


22.60 


21.80 


21.20 


22.70 


22.00 


21.20 


23.90 


23.00 


22.10 


25.30 


24.70 


23.90 


21.10 


21.10 


20.90 


23.60 


23.20 


22.50 



log(Mass) 

[log(M0)] 

(7) 



21.62 
18.51 

20.93 
21.76 








2.88 
1.06 
19.1 
2.14 
3.26 
5.26 


9.52t"oi;^ 
10.05lH« 

9.64«:1«3 

io.25J:S:?| 


20.16 








2.40 


10.22tH? 


20.62 








2.28 


io.o6l«:~ 


20.15 
20.69 








3.15 
1.55 
11.4 


q Qo+0.31 

10.12+«-30 

10.231?-? 


21.83, 


21.35, 


21.20, 


20.59 


4.22 


10.20t«:^^ 


21.93, 


21.55, 


21.50, 


20.89 


3.52 


iom+_l:\l 


24.43, 


23.85, 


23.60, 


21.59 


5.31 


9.93t°:g« 


22.43, 


21.85, 


21.60, 


20.80 


2.01 


io.i2ig:^e 


22.03, 


21.65, 


21.40, 


20.59 


2.65 


io.iil«:io 


21.63, 


21.65, 


21.40, 


21.09 


5.69 


smtVil 


23.93, 


23.45, 


23.20, 


21.59 


1.12 


-.r^^O + O-Ol 
iU.Z/_(J Qg 


25.35, 


24.54, 


24.39 




2.14 


o OC--I-0.37 
o.oo_Q 28 


24.20, 


23.80 






13.0 


8.63«:lt 


23.10, 


23.00 






5.30 


8.89«;1? 


21.40, 


21.10 






22.0 


9.43«-°^9 


23.30, 


23.30 






5.90 


w.i5+_l:lt 


20.20, 


19.80 






8.20 


lo.ostlil 


20.00, 


19.60 






14.0 


9.61«:«« 


20.20, 


19.90 






19.0 


9.59«i| 


21.00, 


20.60 






5.40 


10.19t«:i? 


20.80, 


20.30 






9.80 


lomtlf. 


21.70, 


21.20 






20.0 


10.26t?,:i 


23.00, 


22.90 






13.0 


io.22iEJ:i? 


20.50, 


20.20 






73.0 


q qo+0.04 


22.20, 


21.70 






1.90 


q 74+0.36 



/x//r 
(8) 



1.46 
0.32 
0.33 
1.26 
0.56 
2.54 
0.25 
0.37 
0.18 
0.67 
0.84 
0.37 
0.37 
5.15 
0.31 
0.31 
0.51 
0.65 
4.49 
9.77 
2.45 
1.60 
1.21 
0.25 
0.33 
0.64 
0.31 
0.57 
2.55 
7.49 
2.65 
0.32 



CO 



a 
?= 



Note. — Col.(l): Source name. Col(2): Redshift. Col. (3): RA and DEC of the X-ray target. Col. (4): The distance in arcsec of the optical counterpart from the X-ray target. Col. (5): 
Optical/near-IR photometry in AB system. The photometric bands for each field arc listed in Tablc^ Col. (6): The observed-frame 2-8 KeV X-ray flux. Col. (7): The stellar mass. Col.(S): The 
hard X-ray to i?-band flux ratio. 



